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Abstract 

The e + -jet invariant mass spectrum produced in the reaction e + p — > e + X has been 
studied at a center-of-mass energy of 300 GeV. The data were collected using the 
ZEUS detector operating at the HERA collider, and correspond to an integrated 
luminosity of 47.7 pb _1 . The observed mass spectrum is in good agreement with 
Standard Model expectations up to an e + -jet mass of 210 GeV. Above this mass, 
some excess is seen. The angular distribution of these events is typical of high-Q 2 
neutral current events and does not give convincing evidence for the presence of a 
narrow scalar or vector state. Limits are presented on the product of cross section 
and branching ratio for such a state and are interpreted as limits on leptoquark 
or R-parity-violating squark production. Specific leptoquark types are ruled out 
at 95% confidence level for coupling strength A = 0.3 for masses between 150 and 
280 GeV. 
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1 Introduction 



The e + -jet mass spectrum in e + p scattering has been investigated with the ZEUS detec- 
tor at HERA. An excess of events relative to Standard Model expectations has previously 
been reported by the HI [lj, @ an d ZEUS ||, f§ collaborations in neutral current deep 
inelastic scattering at high x and high Q 2 . These events contain high-mass e + -jet final 
states. Several models have been discussed || as possible sources of these events, includ- 
ing leptoquark production |J and R-parity- violating squark production |7j]. This paper 
presents an analysis of ZEUS data specifically aimed at searching for high mass states 
decaying to e + -jet. 

Candidate events with high transverse energy, an identified final-state positron, and at 
least one jet are selected. The measured energies (E' e , Ej) and angles of the final-state 
positron and the jet with highest transverse momentum are used to calculate an invariant 
mass 

M 2 ej = 2E' e E r (l-cosO, (1) 

where £ is the angle between the positron and jet. The angle between the outgoing and 
incoming positron in the e + -jet rest frame, 8*, is also reconstructed using the measured 
energies and angles. No assumptions about the production process are made in the 
reconstruction of either M e j or 6*. 

The search was performed using 47.7 pb" 1 of data collected in the 1994-1997 running 
periods. In the following, expectations from the Standard Model, leptoquark produc- 
tion and R-parity-violating squark production are summarized. After a discussion of the 
experimental conditions, the analysis is described and the M e j and cos#* distributions 
presented. Since these distributions do not show a clear signal for a narrow resonance, 
limits on the cross section times branching ratio are extracted for the production of such 
a state. Limits are also presented in the mass versus coupling plane which can be applied 
to leptoquark and squark production. 



2 Model Expectations 

High-mass e + -jet pairs, produced in the Standard Model (SM) via neutral current (NC) 
scattering, form the principal background to the search for heavy states. This process is 
reviewed first. Leptoquark (LQ) production and squark production in R-parity-violating 
{ftp) supersymmetry are used as examples of physics beyond the SM that could generate 
the e + -jet final state. The diagrams for NC and LQ processes are shown in Fig. 0. The 
squark production diagrams are similar to the LQ diagrams, but different decay modes 
are possible, as discussed below. 

2.1 Standard Model Expectations 

The kinematic variables used to describe the deep inelastic scattering (DIS) reaction 

e + p — > e + X 
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Figure 1: Diagrams for NC scattering via a) photon exchange and b) Z° exchange. The 
leptoquark diagrams for the same initial and final states are c) s-channel LQ production 
for fermion number F = LQ and d) u-channel LQ exchange for an F = 2 LQ. 



are 

,2 _ fu 1J\2 



Q = -q 2 = -(k-k'Y } (2) 
q-P 



k-P 

Q 2 



and (3) 



1 = Wp' (4 > 

where k and k' are the four-momenta of the incoming and outgoing positron, respectively, 
and P is the four-momentum of the incoming proton. The center-of-mass energy is given 
by s = (k + P) 2 w (300 GeV) 2 . The NC interaction occurs between the positron and a 
parton (quark) inside the proton (see Fig. [!]). The production of the large e + -jet masses of 
interest requires high x partons, where the valence quarks dominate the proton structure. 

In leading-order electroweak theory, the cross section for the NC DIS reaction can be 
expressed as [§] 

^A = ^[ Y+Fi -Y_ xFsWFl \ (5 ) 
dxdy sx A y A L J 

with Y± = 1 ± (1 — y) 2 and a the fine structure constant. The contribution from the 
longitudinal structure function, F L , is expected to be negligible in the kinematic range 
considered here. 

The x dependence of the NC cross section is very steep. In addition to the explicit 1/x 2 
factor, the structure functions F 2 and xF s are dominated at large x by valence-quark 
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densities that fall quickly for x > 0.3. The y dependence of the cross section is dominated 
by the 1/y 2 term. The structure functions vary slowly with y at fixed x. The uncertainty 
in the NC cross section predicted by Eq. || is dominated by the uncertainty in the structure 
functions (parton densities), and is small, about 5% at the high-x and moderate-y ranges 
of this analysis (§]. The quantity of interest in this paper is the e + -jet cross section, which 
is sensitive to QCD corrections. The uncertainty arising from these corrections has been 
estimated to be small for this analysis. 

For DIS or LQ events produced via the diagrams shown in Fig. [I] (i.e. assuming no QED 
or QCD radiation), the mass of the eq system is related to x via 

M 2 = sx (6) 

and 9* is related to y via 



cos0* = l-2y. (7) 

The steeply-falling x and y dependences of DIS events will therefore produce distributions 
falling sharply with mass and peaking towards cos^* = 1. 



2.2 Leptoquark Production and Exchange 

Leptoquark production is an example of new physics that could generate high-mass e + -jet 
pairs. The set of leptoquarks with SU(3) x SU(2) x £/(l)-invariant couplings has been 
specified |5| . Only LQs with fermion number F = L + 3B = are considered here, where 
L and B denote the lepton and baryon number, respectively. These leptoquarks are listed 
in Table |T] together with some of their properties. The F = LQs have higher cross 
sections in e + p scattering than e~p scattering since in the e + p case a valence quark can 
fuse with the positron. 

In principle, additional LQ types can be defined || which depend on the generations of 
the quarks and leptons to which they couple. Only LQs which preserve lepton flavor and 
which couple to first-generation quarks are considered in this analysis. 

As shown in Fig. [1], leptoquark production can generate an s-channel resonance provided 
m LQ < \/s. Contributions to the e + p cross section would also result from u-channel 
exchange and interference of LQ diagrams with photon and Z° exchange. The cross 
section in the presence of a leptoquark can be written as 



dxdy dxdy dxdy dxdy dxdy dxdy 

The first term on the right-hand side of Eq. |8| represents the SM contribution discussed 
previously. The second (third) term arises from the interference between the SM and u- 
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LQ species 


q 


Production 


Decay 


Branching ratio 


qL 

J./ 4 


-5/3 


e L u 


eu 


1 


cR 


-5/3 


e R u 


eu 


1 


-2/3 


e R d 


ed 


1 


1/2 


-2/3 


e L d 


ed 


1 


v L 


-2/3 


e L d 


ed 
v e u 


1/2 
1/2 


v R 


-2/3 


e R d 


ed 


1 


v« 


-5/3 


e R u 


eu 


1 




-5/3 


e L u 


eu 


1 




-2/3 


e L d 


ed 
v e u 


1/2 
1/2 



Table 1: The F = leptoquarks that can be produced at HERA. The LQ species are divided 
according to their spin (S for scalar and V for vector), their chirality (L or R) and their 
weak isospin (0, 1/2,1). The leptoquarks S and V differ by two units of hypercharge from S 
and V, respectively. In addition, the electric charge, q, of the leptoquarks, the production 
channel, as well as their allowed decay channels assuming lepton-flavor conservation, 
are displayed. The quantum numbers and decay channels correspond to an electron-type 
LQ. For positrons, the corresponding anti-leptoquarks have the sign of the electric charge 
reversed, the helicity of the incoming lepton reversed and antiquarks are replaced by the 
corresponding quark. The nomenclature follows the Aachen convention |77|/ . 

channel (s-channel) LQ diagram, and the fourth (fifth) term represents the u-channel (s- 
channel) LQ diagram alone. The additional contributions to the SM cross section depend 
on two parameters: m L Q, and \ R or X L , the coupling to e\ R and quark. Leptoquarks of 
well-defined helicity [\ r -\l = 0) are assumed for simplicity in the limit-setting procedure, 
and one species of LQ is assumed to dominate the cross section. The cos 8* dependence 
varies strongly for the different terms: it is flat for scalar-LQ production in the s-channel, 
and for vector-LQ exchange in the u-channel, while it varies as (1 + cos#*) 2 for vector-LQ 
production in the s-channel or scalar-LQ exchange in the u-channel. The interference 
terms produce a cos 9* dependence which is steeper due to the sharply-peaking cos#* 
distribution in NC DIS. 

In general, the s-channel term dominates the additional contributions to the SM cross 
section if ttilq < \fs, the coupling A is small, and the LQ is produced from a quark rather 
than an antiquark. However, there are conditions for which the other terms can become 
significant, or even dominant Oj, leading to important consequences for the expected 
mass spectra and decay angular distributions. The u-channel and interference terms 
cannot produce a resonance peak in the mass spectrum and the angular distributions 
from such terms can behave more like those of NC deep inelastic scattering. Limits are 
presented in this paper for narrow-width LQ and under conditions for which the s-channel 
term dominates. 
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The width of a LQ depends on its spin and decay modes, and is proportional to m L Q 
times the square of the coupling. In the narrow-width approximation, the LQ production 
cross section is given by integrating the s-channel term ||: 

a NWA = {J+l)^ A 2 g(* ,/i) (9) 

where J represents the spin of the LQ, q(xo, fi) is the quark density evaluated at xq = 
m 2 L Q/s and with the scale \i = m 2 L Q. In the limit-setting procedure (Sect. H), this cross 
section was corrected for expected QED and QCD (for scalar LQ only) radiative effects. 
The QCD corrections |12j enhance the cross section by 20 - 30% for the F = LQ 
considered here. The effect of QED radiation on the LQ production cross section was 
calculated and was found to decrease the cross section by 5-25% as ttilq increases from 
100 -> 290 GeV. 



2.3 R-Parity- Violating Squark Production 

In the supersymmetry (SUSY) superpotential, i?-parity-violating terms of the form 
XijkL % L Q 3 L D R [0 are of particular interest for lepton-hadron collisions. Here, L L , Q L , 
and Dr denote left-handed lepton and quark doublets and the right-handed down-type 
quark-singlet chiral superfields, respectively. The indices i, j, and k label their respective 
generations. 

For i = 1, which is the case for ep collisions, these operators can lead to u- and <i-type 
squark production. There are 9 possible production couplings probed in e + p collisions, 



corresponding to the reactions [13 



e + + Uj — > d k , (10) 
e + + d k — > Uj . (11) 

For production and decay via the A' ljfe coupling, squarks behave like scalar leptoquarks 
and the final state is indistinguishable, event by event, from Standard Model neutral and 
charged current events. However, as for the scalar leptoquarks, the angular distributions 
of the final-state lepton and quark will be different and this fact can be exploited in 
performing searches. Limits derived for scalar LQ production can then be directly related 
to limits on squark production and decay via A- jfc . In addition to the Yukawa couplings, 
gauge couplings also exist whereby q can decay by radiating a neutralino or chargino 
which can subsequently decay. The final-state signature depends on the properties of the 
neutralino or chargino. The search for such decay topologies from a squark is outside the 
scope of this analysis. 



3 Experimental Conditions 

In 1994-97, HERA operated with protons of energy E p = 820 GeV and positrons of 
energy E e = 27.5 GeV. The ZEUS detector is described in detail in |TJ[]. The main com- 
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ponents used in the present analysis were the central tracking detector (CTD) positioned 
in a 1.43 T solenoidal magnetic field and the uranium-scintillator sampling calorimeter 
(CAL). The CTD was used to establish an interaction vertex with a typical resolution 
along (transverse to) the beam direction of 0.4 (0.1) cm. It was also used in the positron- 
finding algorithm that associated a charged track with an energy deposit in the calorime- 
ter. The CAL was used to measure the positron and hadronic energies. The CAL consists 
of a forward part (FCAL), a barrel part (BCAL) and a rear part (RCAL), with depths of 
7, 5 and 4 interaction lengths, respectively. The FCAL and BCAL are segmented longitu- 
dinally into an electromagnetic section (EMC), and two hadronic sections (HAC1,2). The 
RCAL has one EMC and one HAC section. The cell structure is formed by scintillator 
tiles; cell sizes range from 5 x 20 cm 2 (FEMC) to 24.4 x 35.2 cm 2 at the front face of a 
BCAL HAC2 cell. The light generated in the scintillator is collected on both sides of the 
module by wavelength-shifter bars, allowing a coordinate measurement based on knowl- 
edge of the attenuation length in the scintillator. The light is converted into an electronic 
signal by photomultiplier tubes. The cells are arranged into towers consisting of 4 EMC 
cells, a HAC1 cell and a HAC2 cell (in FCAL and BCAL). The transverse dimensions of 
the towers in FCAL are 20 x 20 cm 2 . One tower is absent at the center of the FCAL and 
RCAL to allow space for passage of the beams. The outer boundary of the inner ring 
of FCAL towers, used to define a fiducial cut for the jet reconstruction, defines a box of 
60 x 60 cm 2 . 

Under test beam conditions, the CAL has energy resolutions of a/E = 0.18/ \/E~ for 
positrons hitting the center of a calorimeter cell and a/E = 0.35/v^E for single hadrons, 
where energies are in GeV. In the ZEUS detector, the energy measurement is affected 
by the energy loss in the material between the interaction point and the calorimeter. 
For the events selected in this analysis, the positrons predominantly strike the BCAL, 
while the jets hit the FCAL. The in-situ positron-energy resolution in the BCAL has 
been determined to average a/E = 0.32 /yE © 0.03 while the jet-energy resolution in 
the FCAL averages a/E = 0.55/ \fE © 0.02. The jet-energy resolution was determined 
by comparing reconstructed jet energies in the calorimeter with the total energy of the 
particles in the hadronic final-state using Monte Carlo simulation, and therefore includes 
small contributions from the jet-finding algorithm. 

In the reconstruction of the positron and jet energies, corrections were applied for inactive 
materials located in front of the calorimeter and for non- uniformities in the calorimeter 
response 0]. For the high energies important in this analysis, the overall energy scale is 
known to 1% for positrons in BCAL and 2% for hadrons in FCAL and BCAL. The elec- 
tromagnetic energy scale was determined by a comparison with momentum measurements 
in the central tracking detector (using lower-energy electrons and positrons). Its linearity 
was checked with energies reconstructed from the double angle (DA) method JOJ. The 
hadronic-energy scales in the FCAL and BCAL were determined by using transverse- 
momentum balance in NC DIS events. 

The angular reconstruction was performed using a combination of tracking and calorimeter 
information. From Monte Carlo studies, the polar-angle resolutions were found to be 
2.5 mrad for positrons and approximately (220/ \[E — 4) mrad for jets with energies 
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above 100 GeV. 



The luminosity was measured from the rate of the bremsstrahlung process e + p —>■ e + jrj |16 
and has an uncertainty of 1.6%. 

The ZEUS coordinate system is right-handed and centered on the nominal interaction 
point, with the Z axis pointing in the direction of the proton beam (forward) and the X 
axis pointing horizontally toward the center of HERA. The polar angle 9 is defined with 
respect to the Z axis. 



4 Event Selection 



The events of interest with large e + -jet mass contain a final-state positron at a large angle 
and of much higher energy than that of the incident positron beam, as well as one or more 
energetic jets. The only important SM source of such events is NC scattering with large 
Q 2 . Other potential backgrounds, such as high transverse-energy (E?) photoproduction, 
were determined to be negligible. 

The following requirements selected events of the desired topology: 



A reconstructed event vertex was required in the range \Z\ < 50 cm. 

The total transverse energy, E?, was required to be at least 60 GeV. 

An identified Q positron was required with energy E' e > 25 GeV, located either 
in the FCAL or BCAL. The positron was required to be well-contained in the 
BCAL or FCAL and not to point to the BCAL/FCAL interface, at approximately 
31° < 9 < 36°. Positrons within 1.5 cm of the boundary between adjacent BCAL 
modules, as determined by tracking information, were also discarded to remove 
showers developing in the wavelength-shifter bars. 

A hadronic jet with transverse momentum > 15 GeV, located in a region of good 
containment, was required. The jets were reconstructed using the longitudinally- 



invariant /^-clustering algorithm |17j in the inclusive mode [fig] . Only jets with 
a reconstructed centroid outside the inner ring of FCAL towers were considered. 
In events where multiple jets were reconstructed, the jet with highest transverse 
momentum was used. After all cuts, 12% of the events had more than one jet, both 
in the data and Monte Carlo simulation (see below). 



The Et cut, the jet- containment cut and the positron-containment cut define the available 
kinematic region for further analysis, as shown in Fig. ^. The jet containment cut, in 
particular, limits the values of cos 9* that can be measured at the highest e + -jet masses. 
Because most such events have cos#* near 1, the acceptance for NC DIS events (with 
E T > 60 GeV) falls below 10% for masses beyond 220 GeV. In the region allowed by the 
cuts shown in Fig. |2|, the acceptance is typically 80%. 
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50 100 150 200 250 300 

Mass(GeV) 



Figure 2: The acceptance region (unshaded) in the cos 9* versus M e j plane allowed by the 
Et, jet- containment and positron-fiducial-volume cuts, assuming eq — * eq scattering at 
the nominal interaction point. No detector simulation is included. 
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A total of 7103 events remained after applying all cuts, compared to 6949 ±445 events pre- 
dicted by the NC Monte Carlo simulation based on the measured luminosity of 47.7 pb _1 



(the sources of uncertainty on the expected number of events are described in Sect. |7TT 



The Et distributions for data and NC simulation are compared in Fig. [3]a. The positron 
transverse-momentum (P T ) spectrum, jet transverse-momentum (P T ) spectrum and the 
ratio P^/Pti where P? is the transverse momentum of the hadronic system, are shown in 
Figs, ^(b-d), respectively. The missing transverse momentum, P?, and the longitudinal 
momentum variable, E — Pz , are compared in Figs. |||(e,f). The global properties of the 
events are well reproduced by the simulation. 



5 Event Simulation 



The SM deep inelastic scattering events were simulated using the HERACLES 4.5.2 [19 



program with the DJANGO 6 version 2.4 |20| interface to the hadronization programs. 
In HERACLES, corrections for initial- and final-state electroweak radiation, vertex and 
propagator corrections, and two-boson exchange are included. The NC DIS hadronic final 
state was simulated using the MEPS model of LEPTO 6.5 plj, which includes order as 



matrix elements and models of higher-order QCD radiation. As a systematic check, the 
NC final state was simulated using the color-dipole model of ARIADNE 4.08 |[22|| . The 
CTEQ4 parton-distribution set |23) was used to evaluate the expected number of events 
from NC DIS scattering. 



The leptoquark events were generated using PYTHIA 6.1 [21]. This program takes into 
account the finite width of the LQ, but only includes the s-channel diagrams. Initial- and 
final-state QCD radiation from the quark and the effect of LQ hadronization before decay 
are taken into account, as are initial- and final-state QED radiation from the positron. 



The generated events were input into a GEANT-based |25j program which simulated the 
response of the ZEUS detector. The trigger and offline processing requirements applied 
to the data were applied to the simulated events. The luminosity of the NC Monte Carlo 
samples ranges from 46 pb" 1 at Q 2 = 400 GeV 2 to 7.3 ■ 10 6 pb~ x at Q 2 = 50000 GeV 2 . 



6 Mass and 6* Reconstruction 

The mass of each e + -jet pair was reconstructed from the measured energies and angles of 
the positron and jet as described by Eq. |I[ This formula makes no correction for the finite 
jet mass. Possible mass shifts and the resolutions for resonant lepton-hadron states were 
estimated from PYTHIA. Narrow scalar LQ events in the mass range 150 — 290 GeV were 
simulated. The mean mass for reconstructed events was found to be within 6% of the 
generated value, while the peak position as determined by a Gaussian fit was typically 
lower than the generated value by only 1%. The average mass resolution, determined 
from a Gaussian fit to the peak of the reconstructed mass spectrum, ranged from 5.5% to 
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Figure 3: Comparison of data (points) with Standard Model expectations (histograms) for 
selected distributions: a) total transverse energy, Et; b) positron transverse momentum, 
Pt'j c ) ) e t transverse momentum, Pj>; d) the ratio of the positron to hadron transverse 
momenta, P^/P?, (e) the missing transverse momentum, Pt, and (f) the longitudinal- 
momentum variable, E — P z , for the event. 
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3% for masses from 150 to 290 GeV. The RMS of the distribution was typically twice as 
large. 

The positron scattering angle in the e + -jet rest frame, 9*, was reconstructed as the angle 
between the incoming and outgoing positron directions in this frame. These directions 
were determined by performing a Lorentz transformation using the measured positron and 
jet energies and angles in the laboratory frame. The resolution in cos 9* near | cos^*! = 1, 
as determined from a Gaussian fit, was 0.01 degrading to 0.03 as | cos#*| decreases. The 
shift in cos9* was less than 0.01 for both the NC MC and the leptoquark MC. 

In order to determine limits on leptoquark and squark production, the mass of the 
electron-hadron system was reconstructed by the constrained-mass method. This method 
reconstructs the e + -hadron mass as 



where (E + Pz) is the sum of the energy and Pz contributions from the positron and 
all jets satisfying Pf > 15 GeV and pseudorapidity rjj < 3 (with the highest jet 
required to be outside the FCAL inner ring). The rjj cut removes contributions from 
the proton remnant. The constraints fx = and E — Pz = 2P e , which are satisfied by 
fully contained events, have been assumed in arriving at this equation. When using this 
mass-reconstruction method, events with measured E — Pz < 40 GeV were removed to 
avoid large initial-state QED radiation. 

The Mcj method gave, on average, improved resolution over the M e j method for narrow 
LQ MC events. The improved resolution occurred at smaller cos#* (for cos 9* ~ the 
mass resolution determined from a Gaussian fit to the reconstructed mass distribution 
for rriiQ = 200 GeV was about 1.5% in the Mcj method and 3% for the M ej - method); 
at the larger cos#* values where NC DIS events are concentrated, the resolutions of the 
two methods were similar (about 3% for ttilq = 200 GeV). The Mqj method relies on 
constraints which do not necessarily apply to a resonant state whose properties cannot be 
anticipated in detail. We therefore choose to use the M e j method as our primary search 
method. The Mqj method is used in the limit-setting procedure. 

7 M e j and cos#* Distributions 

The reconstructed values of M e j are plotted versus cos 9* for the selected events in Fig. £|. 
Most of the events are concentrated at large cos#* and small mass, as expected from 
Standard Model NC scattering. The five events indicated as open circles are from data 
taken in 1994-96, with total luminosity 20 pb _1 . They were the subject of a previous 
publication ||. In this earlier analysis, the kinematic variables were reconstructed with 
the DA method. The five events also stand out with the M e j reconstruction technique. 
The average value of M e j for these events is 224 GeV, or 7 GeV less than the corresponding 
mass calculated previously via M = ■ Xda, where xda is the estimator of Bjorken-x 
calculated with the DA method. This mass shift is compatible with expectations based 




(12) 



11 



<x> 



o 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 







ZEUS 1994-97 




50 



100 



150 



200 



250 



300 



M ej (GeV) 



gure 4: The value of cos 9* versus M e j for all 7103 events passing the event selections, 
e events at high x and y described in a previous paper J^j are shown as open circles. 
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on resolution and initial state radiation effects. With the present luminosity of 47.7 pb _1 , 
7 events are observed in the region of M e j > 200 GeV and cosO* < 0.5, where 5.0 events 
are expected. 

The M e j spectrum for events with M e j > 100 GeV is shown in Fig. ^a on a logarithmic 
scale. The high-mass part of the spectrum is shown on a linear scale in the inset. The 
predicted number of events ( y N pred ) from NC processes is shown as the histogram. The 
ratio of the measured mass spectrum to the expectation is shown in Fig. |5|b. The shaded 
band indicates the systematic uncertainty on the expectations. 
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Figure 5: a) Comparison of observed events (points) and SM expectations (histogram) for 
the reconstructed e + -jet invariant mass. The inset shows the region with M e j > 180 GeV 
on a linear scale, b) The ratio of the number of observed events to the Standard Model 
expectations. The shaded band shows the systematic uncertainty in the predicted number 
of events. The error bars on the data points are calculated from the square root of the 
number of events in the bin. 
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7.1 Systematic Uncertainties 



The uncertainty on N pred varies with mass from 7% at 100 GeV up to 30% at 250 GeV. 
The most important uncertainties are on the energy scale and the jet position. The NC 
DIS cross section given in Eq. |5| (neglecting F^) can be rewritten in terms of the e + q 
invariant mass, M, and the polar angle of the outgoing struck quark in the laboratory 
frame, 7: 

dV(e+p) 32na^ e sin 7 

~dMd^ ~ MHl ~ cos 7 ) 2 [ + 2 " 3] ' (13) 

The mass dependence is very steep. In addition to the explicit M~ 5 dependence, there is 
also a strong suppression of high masses implicit in the structure functions. An incorrect 
energy scale will produce a shift in the mass spectrum and potentially a significant error 
in the number of expected events at a given mass. The dependence on the quark angle is 
also steep, approximately 7~ 3 at small 7. The number of events passing the jet fiducial 
cut is therefore strongly dependent on the accuracy of the jet position reconstruction. 
The jet fiducial-volume cut requires the highest-Py jet to point outside the inner ring 
of FCAL towers. Many distributions from data and MC were compared to search for 
possible systematic biases. 

The dominant sources of uncertainty are itemized below in order of decreasing importance: 



1. Knowledge of the calorimeter energy scales: 

The scale uncertainties discussed in Sect. |3] are 1% for BCAL positrons and 2% for 
hadrons, leading to an uncertainty of 5(18)% in NP red at M ej = 100(210) GeV. 

2. Uncertainties in the simulation of the hadronic energy flow, including simulation of 
the proton remnant, the energy flow between the struck quark and proton remnant, 
and possible detector effects in the innermost calorimeter towers: 

Many distributions of data and MC were compared and no important systematic 
differences were found. Figure |6| shows the fraction of the jet energy in the inner 
ring of FCAL towers associated with the highest Pt jet as a function of rjj. This is 
shown for all events in Fig. ||a, as well as for those with M e j > 210 GeV in Fig. |^b. 
For the highest rjj values considered, this ratio is about 20%. The energy located 
in the innermost towers of the FCAL and not associated with the highest Pt jet 
is shown in Fig. ||c,d, and compared to the MC simulation. No large differences 
are seen between data and MC (the lowest 77 bin in Fig |]d contains only five data 
events). The innermost towers of the FCAL have a larger uncertainty in the energy 
scale than the rest of the FCAL owing to their slightly different construction and 
proximity to the beam. The energy in these cells has been varied by ±10%. As 
a test of the simulation of the forward energy flow, the ARIADNE MC has been 
used instead of the LEPTO MC. These tests yielded variations in N pred of 13% at 
M ej = 210 GeV. 

3. Uncertainty in the parton density functions: 
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Figure 6: The ratio of the jet energy in the innermost towers of the FCAL, EIRj, to 
the total jet energy, Ej, as a function of r] of the jet for a) the full sample, and b) for 
those events with M e j > 210 GeV. The energy deposited in the innermost FCAL towers, 
excluding that associated with the highest jet, is shown in c) for the full sample, and 
in d) for those events with M e j > 210 GeV. The data are shown as points, while the NC 
Monte Carlo predictions are shown as a histogram. 
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The parton density functions were estimated as in Q, and led to an uncertainty of 
5% in NP red at M ej = 210 GeV. 

4. Uncertainties in the acceptance: 

The alignment of the FCAL was determined to better than 5 mm, and various 
jet position reconstruction algorithms were compared. These studies yielded an 
uncertainty of 2% in N pred . 

5. Uncertainties in the energy resolution functions. 

These were studied by comparing tracking information with calorimeter information 
for individual events, as well as by comparing different reconstruction methods. The 
MC energies were smeared by additional amounts to represent these uncertainties, 
leading to a variation of less than 5% in N pred . 

Other uncertainties include positron finding efficiency, luminosity determination, vertex 
simulation, multijet production rates, and hadronization simulation. These were found to 
be small in comparison to the items listed above. The overall systematic uncertainty was 
obtained by summing the contributions from all these sources in quadrature. 

7.2 Discussion 

The data in Fig. [| are in good agreement with the SM expectations up to M e j « 210 GeV. 
Some excess is seen at higher masses. For M e j > 210 GeV, 49 events were observed in 
the data, while 24.7 ± 5.6 events are expected. A careful study of individual events in 
this mass region uncovered no signs of reconstruction errors. Rather, the events always 
contain clear examples of a high-energy positron (typically 70 GeV) near 90° and a high- 
energy jet (typically 400 GeV) in the forward direction (2 events have a second jet, in 
accord with NC DIS Monte Carlo expectations). The distributions shown in Fig. § for 
all the data are restricted to the events with M e j > 210 GeV in Fig. [7]. Whereas the 
shapes of the distributions are similar, the data lie systematically above the MC, which 
is normalized to the integrated luminosity. 

The events with large M e j have characteristics similar on average to NC DIS events. In 
particular, the cos 8* projection of the events with M e j > 210 GeV is shown in Fig. |8| 
and compared to the MC expectations for neutral current DIS (solid histogram). The 
expectations for narrow s-channel scalar and vector LQ production are also shown for 
comparison. For F = LQs with A < 1, the u-channel and interference terms would 
not significantly affect these expectations. The shape of the data and NC MC cos#* 
distributions are qualitatively similar, peaking at high values of cos 8*. 

In summary, there is some excess of events with M e j > 210 GeV above the Standard 
Model predictions. The probability of observing such an excess depends strongly on 
possible systematic biases. The most important of these are biases in the energy scales. 
As a test, many MC experiments were generated where the jet energy scale was shifted by 
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Figure 7: Comparison of data (points) with Standard Model expectations (histograms) for 
selected distributions and requiring M e j > 210 GeV: a) total event transverse energy, Et! 
b) positron transverse momentum, P^; c) jet transverse momentum, P?>; d) the ratio of 
the positron to hadron transverse momenta, P^/P^, (e) the net (or missing) transverse 
momentum, Pt, and (f) the longitudinal momentum variable, E — P z , for the event. 
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+2% and the electron energy scale by +1%. A window of width 3a(M e j), where a(M e j) 
is the mass resolution at mass M ej -, was moved over the accessible mass range. For 
each simulated experiment, the number of observed events within the mass window was 
compared with the nominal expectations as a function of M e j, seeking the excess which 
gave the largest statistical significance. The same procedure was applied to the data. As 
a result, it was found that 5% of the simulated experiments would observe, somewhere 
in the mass spectrum, an excess of statistical significance at least as large as the one 
found in the data. The excess is therefore not statistically compelling. Furthermore, the 
events have the characteristics of neutral current scattering. Limits are therefore set on 
the production of narrow scalar or vector states, as discussed below. 
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Figure 8: Comparison of data (points) and NC Monte Carlo expectations (solid histogram) 
for the cos 8* distribution for events with M e j > 210 GeV. The predictions for narrow 
scalar and vector resonant leptoquarks are shown with arbitrary normalization for com- 
parison. 



8 Limits on Narrow Scalar and Vector States 

Limits are set on the production cross section times branching ratio into positron+jet, 
aB, for a narrow scalar or vector state. For definiteness, limits on coupling strength 
versus mass for F = leptoquarks are presented, as well as limits on X^/B versus mass 
for scalar states coupling to u or d quarks, such as ftp squarks. The limits are extracted 
for A < 1, allowing the use of the narrow- width approximation assumed in Eq. |9|. 

The Mc.j mass reconstruction method was used to set limits as described in Sect. [5[ The 
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Figure 9: The reconstructed mass spectrum using the Mqj method for data (points) and 
SM expectations (histogram) : a) shows the spectrum for events satisfying all cuts, while 
b) shows the mass spectrum after the cut (cos 6* < cos$* ut ) for the scalar LQ search has 
been applied. 

positron fiducial cuts were removed since this method is less sensitive to the positron- 
energy measurement, while the cut E — Pz > 40 GeV was applied to reduce radiative 
effects. The mass spectrum reconstructed with this technique is shown in Fig. |)]a. In 
total, 8026 events passed all selection cuts while 7863 events are predicted by the MC. 

The leptoquark MC described in Sect. [5] was used to determine the event selection effi- 
ciency and the acceptance of the fiducial cuts, as well as to estimate the mass resolution. 
This MC and the NC background simulation were used to calculate an optimal bin width, 
AM C j, for each M C j, and optimal cos^* range, cos6>* < cos9* ut , to obtain on average 
the best limits on LQ couplings. The bin widths were typically 20 GeV. The values of 
cos6* ut for setting limits range from 0.5 to 0.9 for vector leptoquarks with masses between 
150 — 290 GeV, and from 0.1 to 0.9 for scalar leptoquarks in the same mass range. The 
mass spectrum after applying the optimal cos#* cut for the scalar search is shown in 
Fig. ^jo. No significant deviations from expectations are seen after applying this cut. 

The 95% confidence level (CL) limits on o~B were obtained directly from the observed 
number of data events with cos#* < cos8* ut in the particular mass window p6 . The 



procedure described in [^6| was extended to include the systematic uncertainties in the 
numbers of predicted events. This was found to have negligible effect on the limits. The 
limits for a narrow scalar or vector state are shown in Fig. These limits lie between 1 
and 0.1 pb as the mass increases from 150 to 290 GeV. 
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Figure 10: Limits on the production cross section times branching ratio for decay into 
e + -jet(s) for a scalar or vector state, as a function of the mass of the state. The shaded 
regions are excluded. 
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The 95% CL exclusion limits for different species of LQ are given in the coupling versus 
mass plane in Fig. [11]. The limits exclude leptoquarks with coupling strength A = ^/Ana ~ 
0.3 for masses up to 280 GeV for specific types of F=0 leptoquarks. The HI collaboration 
has recently published similar limits D. In Fig. D], the ZEUS results are compared to 
recent limits from OPAL. At LEP [^7- 29f] , sensitivity to a high-mass LQ arises from effects 
of virtual LQ exchange on the hadronic cross section. The HERA and LEP limits are 
complementary to Tevatron limits p0|,pl|, which are independent of the coupling A^. 



The limits by DO (CDF) extend up to 225 (213) GeV for a scalar LQ with 100% branching 
ratio to eg. The DO limits are shown as vertical lines in Fig. [11]. The Tevatron limits for 
vector LQs are model dependent |32|, but are expected to be considerably higher than for 
scalar LQs. 

The ZEUS limits presented in Fig. [11] can also be applied to any narrow state which couples 
to a positron and aaorrf quark and with unknown branching ratio to e + -jet(s). These 
states correspond to the leptoquark types as labelled in the figure. For these states, the 
limits are on the quantity \\fB. Examples of scalar states for which these limits apply 
are /£p-squarks (e.g., the limit on the Sy 2 (e + d) LQ can be read as a limit on the A^- x 
R-parity- violating coupling). 



9 Conclusion 

Data from 47.7 pb _1 of e + p collisions at a center-of-mass energy of 300 GeV have been 
used to search for a resonance decaying into e + -jet. The invariant mass of the e + -jet pair 
was calculated directly from the measured energies and angles of the positron and jet. 
This approach makes no assumptions about the production mechanism of such a state. 

The observed mass spectrum is in good agreement with Standard Model expectations up 
to e + -jet masses of about 210 GeV. Above this mass, some excess is seen. The angular 
distribution of these events is typical of high-Q 2 neutral current events and does not give 
convincing evidence for the presence of a narrow scalar or vector state. By applying 
restrictions on the decay angle to optimize sensitivity to a narrow state in the presence 
of NC background, limits have been derived on the cross section times decay branching 
fraction for a scalar or vector state decaying into positron and jet(s). These limits can be 
interpreted, for example, as limits on leptoquark or R-parity-violating squark production. 
Limits on the production of leptoquarks and squarks are presented in the coupling strength 
versus mass plane. At a coupling strength A = 0.3, new states are ruled out at 95% 
confidence level for masses between 150 and 280 GeV. 
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Figure 11: Coupling limits as a function of leptoquark mass for F = leptoquarks. The 
results from this analysis are compared to representative limits from LEP J^/ and the 
Tevatron fi3l\] . The areas above the ZEUS and LEP curves are excluded, while the area to 
the left of the Tevatron line is excluded for scalar leptoquarks with the indicated branching 
ratio to e+jet. The limits on scalars are shown in a) while the limits on vectors are shown 
in b). 
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